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ABSTRACT 

The analysis of stellar abundances for odd-Z Fe-peak elements requires accurate NLTE mod- 
elling of spectral lines fully taking into account the hyperfine structure splitting (HFS) of 
lines. Here, we investigate the statistical equilibrium of Co in the atmospheres of cool stars, 
and the influence of NLTE and HFS on the formation of Co lines and abundances. Significant 
departures from LTE level populations are found for Co i, also number densities of excited 
states in Co ii differ from LTE at low metallicity. The NLTE level populations are used to 
determine the abundance of Co in solar photosphere, log e = 4.95 + 0.04 dex, which is in 
agreement with that in C i meteorites within the combined uncertainties. The spectral lines of 
Co I were calculated using the results of recent measurements of hyperfine interaction con- 
stants by UV Fourier transform spectrometry. For Co ii, the first laboratory measurements of 
hyperfine structure splitting A and B factors were performed. These highly accurate A factor 
measurements (errors of the order of 3-7%) allow, for the first time, reliable modelling of 
Co II lines in the solar and stellar spectra and, thus, a test of the Co i/Co ii ionization equi- 
librium in stellar atmospheres. A differential abundance analysis of Co is carried out for 18 
stars in the metallicity range -3.12 < [Fe/H] < 0. The abundances are derived by method of 
spectrum synthesis. At low [Fe/H], NLTE abundance corrections for Co i lines are as large as 
-1-0.6 . . . + 0.8 dex. Thus, LTE abundances of Co in metal-poor stars are severely underesti- 
mated. The stellar NLTE abundances determined from the single UV line of Co ii are lower by 
~ 0.5 - 0.6 dex. The discrepancy might be attributed to possible blends that have not been ac- 
counted for in the solar Co ii line and its erroneous oscillator strength. The increasing [Co/Fe] 
trend in metal-poor stars, as calculated from the Co i lines under NLTE, can be explained if 
Co is overproduced relative to Fe in massive stars. The models of galactic chemical evolution 
are wholly inadequate to describe this trend suggesting that the problem is in SN yields. 
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1 INTRODUCTION 

Cobalt is one of the most intriguing elements in studies of galactic 
chemical evolution. There is neither an agreement on the overall 
abundance trend of Co in the halo and thin/thick disk, nor a gener- 
ally accepted scenario for the nucleosynthetic production of cobalt. 

Several investigations of Co abundances have been done in 
various metallicity regimes and stellar populations. All of these 
abundance investigations were based on Co i lines and local ther- 
modynamic equilibrium (LTE) for the line formation. According to 
the results o f lOratton & SnedeiJh99lh . Co follows the depletion of 
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Fe dow n to [Fe/H] » -2.5. At lower metallicities. lMcWilliam et"al] 
d 19951) observe a stro ng increase of cobalt a bundance with respect 
to iron by ~ 0.5 dex. iDel Peloso et al.l ( 12005) obtained a different 
trend for [Co/Fe], characterized by a steady increase towards su- 
persolar ratios of ~ +0.2 from [Fe/H] = to -0.8. An excess of 
co balt with respect t o iron at -4 < [Fe /H] < -2.5 was reporte d 



by iLai et alj ll2008l) ICayrel et alj ( |2004|) . and iRvan et al.l ( Il99lh . 
iReddv et all feool l 20061) find supersolar and mildly subsolar Co 
abundances for the thick and thin disk stars, respectively. 

The increasing [Co/Fe] ratio with decreasing metallicity was 
attributed to the Co production in both types of supernovae (SNe) 
but with uncertain relative contributions ( Del Peloso et al. 200J). 
However, this is hard to understand from the point of view of stellar 
nucleosynthesis. Cobalt is an odd-Z element of the Fe-group and is 
produced in a-rich freeze-out in SNe. Thus, its abundance must be 
sensitive to the explosion entropy and neutron excess available to 
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the progenitor star in the explosive Si-burning phase. The neutron 
enrichment is determined by the initial metal content of a star and 
by the previous hydrostatic burning stages. Hence, in metal-poor 
environments we expect to see an odd-even effect characterized by 
a general suppression of the stable odd-Z nuc lei relative to their 
stable even-Z neighbours teurbidge et"ai]|l957h . Still, this effect is 
in conflict with the increasing [Co/Fe] trend in metal-poor stars. 

Thus far, any attempt to describe the abundances of Co 
within the framework of explosive nucleosynthesis in SNe failed. 
Models of galactic chemical evolution (GCE), which use stan- 
dard metallicity-dependent prescriptions for Co production in SNe 
II Iw oosley & Weaver 1995) and SNe la jNomoto et all Il984 
llwamoto etal.l 119991) predict a decreasing or nearly flat [Co/Fe] 
trend jTimmes et all 1 19951 : ISamlanj|l998l : lOoswami & Prantzoj 
120001) . Models with metallicity-independent SN yields give only 
a qualitatively similar behaviour, however large discrepancies 
with spectroscopic abundan ces are seen eit her for the thick dis k 
jGoswami & Prantzosll2000h or for the halo jprancois et al.l2004l) . 

There are several explanations for the disagreement between 
observed Co trends and predictions of chemical evolution models. 
Deficiencies of the latter are hidden among others in the assumed 
shape of the initial mass function (IMF). As shown by. Wang & Silk 
( Il993h , the shape of the IMF introduces a factor of 2 in uncertainty 
to the absolute yield value of an elemen t. The physics of super - 
nova explosions is still poorly understood jThielemann et al.l2007b . 
hence theoretical stellar yields for Co may be erroneous. It is also 
important to note that the majority of GCE models assume a single 
gas phase for the galactic ISM, and o nly a few consider multiphase 
ISM (e.g. lSchonrich & Binnevl2009l) that is crucial for a physically 
accurate description of the galactic chemo-dynamical evolution. 

On the other hand, a close examination of the 'observed' Co 
abundances in the Sun and metal-poor stars reveals that all of 
them rely on several coarse assumptions: LTE in the line forma- 
tion, static ID atmospheric models, and, often, neglect of hyper- 
fine structure splitting (HFS) of lines. The errors in abundances 
caused by these approximations are hard to predict, especially be- 
cause effec ts of NLTE and temperature inhomogeneities are not in- 
dependent ( lAsplundl200"3) . Iron is the only representative of the Fe- 
group for which calculations of 1.5D NLTE line formation with 3D 
convective model atmospheres have been performed. NLTE anal- 
yses with ID model a tmospheres exist for a number of elements, 
specifically for Fe ( see lKom et a l.*2003l a nd references ther ein), Sc 
jZhang et alj|2008l) . and Mn feergeman n & Gehrenll2008h . These 
and other studies prove that the NLTE abundance corrections can 
be significant at low gravities and/or metallicities. Due to the ove- 
rionization, NLTE abundances of Mn in metal-poor stars are larger 
by +0.1 . . . + 0.6 dex than those computed under LTE. 

Knowledge of hyperfine structure (HFS) in odd-Z iron-group 
elements is essential in calculations of line formation and abun- 
dances. The reason is that HFS is not merely a line broaden- 
ing mechanism; it modifies absorption over the entire line pro- 
file, especially for strong and saturated lines. Thus, neglect of 
HFS results in wavelength shifts, incorrectly determined damp- 
ing paramete rs, and erroneous abundances (this last issue is 
addressed in iHolweger & Oertell [l97ll : IProchaska & McWilUaml 
I2OOOI : iDel Peloso et alJl20o'5ir In a strictly differential analysis of 
stars with respect to the Sun, which is usually preferred in order 
to minimize other deficiencies of the stellar atmosphere modelling, 
neglect of HFS for the reference solar spectral lines may lead to 
systematically underestimated abundances in metal-poor stars. The 
accuracy of the adopted HFS factors is also important, as we will 
show in this paper. At present, high resolution measurements of 



HFS over a wide spectral range are possible with high resolution 
Fourier transform spectrometry (FTS). For example, an analysis of 
Co I was undertaken using over one thousand line profiles acquired 
with the FT spectrometer at Imperial College, and yielded HFS A 
splitti ng factors for 297 energy levels, almost all known Co i energy 
levels l |Pickering|ll996h . 

In this paper, we investigate the excitation-ionization equilib- 
rium of Co in stellar atmospheres and predict what errors in abun- 
dance calculations occur when NLTE effects and hyperfine splitting 
of lines are neglected. Given the particular importance of HFS in 
this context, we have measured hyperfine splitting of Co 11 energy 
levels using the technique of FTS. We analyze the abundance of Co 
in the Sun, and compare it with that of meteorites to ascertain the 
reason for the discrepancies. Furthermore, we compute the NLTE 
abundances for a sample of 18 stars, which belong to the thin and 
thick disk, and the halo. The evolution of abundance ratios in the 
Galaxy and some implications for the Fe-peak element nucleosyn- 
thesis are discussed. 



2 MEASUREMENTS OF THE HYPERFINE STRUCTURE 
OF CO II LEVELS 

2.1 Simple Theory 

For atoms with non zero nuclear spin / the fine structure levels un- 
dergo splitting because of the hyperfine interaction between the nu- 
cleus and the electrons. In the absence of perturbations the energy 
of the hyperfine structure multiplets is given (iKopfermana ,195 
by: 

2 2/(2/ - 1)7(27 - 1) 

where Wj is the energy of the fine structure level of quantum num- 
ber J, A and B are the magnetic dipole and the electric quadrupole 
hyperfine interaction constants respectively, and K is defined as 

K = F(F +1)- J(J + 1) - /(/ + 1) 

in which F is the quantum number associated with the total angular 
momentum of the atom, 



/ + /: / + /-!; 



\I-J\ 



The selection rules which govern the hyperfine transitions are: 
AF = 0; ±1 but not F = F = 

In addition there are intensity rules ( lKuhn|[l963) : within a hy- 
perfine multiplet the ratio of the sums of the intensities of all tran- 
sitions from two states with quantum numbers F and F' are in the 
ratio of their statistical weights (2F + 1):(2F' + 1). 

This simple theory was used in the fitting of HFS of the Co 11 
transitions in this work. For Co, nuclear spin 1 = 7/2. No perturba- 
tions were evident. 



2.2 Experimental Details 

The spectra used in this work were obtained as part of an exten- 
sive measurement of the Co 11 spectrum, by FTS with cobalt-neon 
and cobalt-argon hollow cathode lamps as sources, in the wave- 
length region 1420 - 33000 A (70422 - 3000 cm"'), resulting in 
a doub ling of the number of classified lines, o ver 200 new energy 
levels jPickering et al.lll998l: |pickering"l99yB), and new calcula- 
tions of Co II oscillator strengths (Raassen et al.lll998h . These pa- 
pers give identified line lists, including the centre-of-gravity (COG) 



Cobalt lines in spectra of the sun and metal-poor stars 3 



wavenumber and intensity of each transition, together with ta- 
bles of energy levels. Most of the Co ii lines exhibit broadening 
due to hyperfine structure. The spectra used for this HFS study 
were acquired in the visible-UV region with the vis-VUV Fourier 
transform (FT) spectrometer at Imperial College, with resolution 
0.037 cm-' (0.006 A) in the visible, 0.05 cm"' (0.003 A) in the 
UV, the individual line components being limited by their Doppler 
widths alone. The light source was a water-cooled hollow cathode 
lamp run in either Ar at 1 mbar or Ne at 2 - 3 mbar, with currents 
of 600 - 700 mA, using a pure cobalt cathode. 

2.3 Hyperfine Structure Analysis 

The determination of the A and B factors for the Co ii energy lev- 
els presented in this paper was carried out using a H FS analysis 
subro utine in the spectrum analysis program Xgremlin (Nave et al. 
Il997i) . The input information comprises the /-value, the /-values of 
both levels, t he known intensity ratios of the transition components 
( lKuhn|[l962h . the widths of the components (having Gaussian and 
Lorenzian contributions), the COG of the pattern and, where possi- 
ble (see below) an A factor for one of the two levels. A large scale 
study of Co II HFS is ongoing at Imperial College, and here we re- 
port a detailed investigation of A factors for the 3c/^('*P)4j a ^Pj ^ j 
levels and for three of the four 3rf^("*F)4p ^'043 1 levels, with im- 
mediate application in this work. 

It was not possible to determine the electric quadrupole hy- 
perfine interaction constant B to great accuracy, because the con- 
tribution of the last term in Eq. 1 was so small, but B factors are 
estimated. The analysis was particularly challenging because no 
previous measurements existed, but with a comprehensive study 
we were able to set a starting point for the fits of line profiles in 
cases where an initial guess at the A factors could be made from a 
characteristic line profile pattern. For example, when one of the A 
factors of the levels involved in a transition is an order of magnitude 
greater than the other A factor, th e line profile ha s a characteristic 
"flag" pattern (examples shown in lPickerind 19961) . These A factors 
could then be used to find A factors of other levels using different 
transitions, building up the number of known A factors. 

In general, for each level studied the HFS A factor was deter- 
mined from the analysis of line profiles for between 3 and 6 differ- 
ent transitions, however in the case of the even levels reported here 
between 7 and 13 line profiles were fitted. 

The resolution of the FT spectrometer is such that the line 
widths are limited solely by the Doppler widths and so, depending 
on the size of the A factors, the overall pattern is usually well re- 
solved into individual components in the IR region, but only partly 
so in the UV. However, even where a set of line profiles did not ap- 
pear to be resolved into individual components, the fitting still gave 
consistent sets of A factors. 

The high number of fitted lines, the fact that the A factors were 
relatively large in magnitude compared to those of higher lying odd 
levels, and good signal-to-noise ratio of studied lines gave reliable 
A factors for the a ^Pj , , levels with uncertainties ranging between 
3 - 7 %. However, fitting of line profiles to find the A factors for 
the four odd z'D° levels was less reliable because of the smaller 
value of the A factor, resulting in almost no visible structure in the 
line profile. A factors for these odd levels should be considered to 
be estimates. 

The values determined for the A and B splitting factors are 
listed in Table[T] and they are the weighted average values from the 
line fits with weighting dependent on th e signal-to-noise ratio of the 
transitions of the individual line fits jSlackwell- Whitehead et all 




38238 38238 38239 38240 38240 

Wavenumber (cm ) 

Figure 1. An example of an observed (heavier trace in plot) and fit- 
ted (dashed trace) line profile, for a Co ii transition recorded by FTS: 
(38238.965 cm"') 2614.3530 A transition, 3rf'C'P)4i a^Pj - 3£(^(''P)4i 



Table 1. Hyperfine structure interaction constants A and B for selected en- 
ergy levels of Co n. N denotes the number of profiles fitted. 
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49.9 
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9.0 


2 


10 
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47039.102 
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47537.362' 














1 


47848.778^ 


<10 











Note: transitions could not be reliably fitted as part of this work to date, 

^ transitions involving z^D° were observed to be narrow, consistent with 
A < 10 mK, but fitting was unreliable, so an upper limit of the A factor 
alone is given here as a guide. 



I2OO5I) . The uncertainties listed in Table [T] are estimated from the 
range of values of each particular A factor, from the standard devi- 
ation. In Table [T] the first four columns give, for each level investi- 
gated, th e configuration, t erm d esignation, /-value and energy re- 
ported in [Pickering et alj ( Il998h . The fifth and sixth columns give 
the A factor and its uncertainty in mK (lmK=0.001cm ') deter- 
mined in this work. The next column lists, where appropriate, ap- 
proximate values of B factors and their uncertainty. The final col- 
umn gives the number of line profiles fitted in the process of finding 
that A factor. An example of a fitted line profile used in this work 
is shown in Fig.[T] 



3 NLTE LINE FORMATION 

3.1 The methods, model atmospheres, and stellar parameters 

The statistical equilibriu m (SE) calculations for C o are carried out 
with the DETAIL code jSutler & Giddingslll985l) , with the radia- 
tive transfer based on the method of accelerated lambda iteration. 
Spectrum synthesis is performed using the SIU code written 
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by J. Reetz. The profiles of Co i and Co ii lines are computed 
with NLTE atomic level populations from SE calculations. Radi- 
ation damping is derived from the oscillator strengths of transi- 
tions. Damping due to quadrat ic Stark and va n der Waals effects fol- 
lows the expressions given bv lUnsoldfigsi) . However, the atomic 
interaction parameter Cs is calculated using the line half-widths 
due to elastic col l isions with H i computed from the theory of 
lAnstee & 0'MaraH l995^. The log Cf, values and other parameteres 
for each Co line investigated in this work are given in Table|2] 

The hyperfine structure of lines is taken into account by su- 
perposing individually synthesised HPS components. Relative in- 
te nsities of the component s are calculated according to the tables 
of 'White &Eliasonl jl933) . Wavelengths are computed using the 
laboratory dat a for A an d B fa ctors. HFS data for all Co i levels 
are taken from lPickerind ( Il996l) . For the levels of Co ii, we use our 
new values as described in Sect. [2] The A and B factors are given 
in the Table lAll available online. 

Abundances are determined by fitting the synthetic lines to 
the observed line profiles. The analysis of stars is strictly differen- 
tia^ relative to the Sun, i.e. any abundance estimate derived from 
a single line in a spectrum of a metal-poor star is referred to that 
from the corresponding solar line. This excludes the use of oscilla- 
tor strengths in calculation of stellar [Co/Fe] ratios. 

All computations, including determinations of stellar param- 
eters, are based on the model atmospheres computed with the 
MAFAGS code. These are static ID plane-parallel models with 
line blanketin g computed with opacity distribution functions of 
lKurucj ( fT993) . Convection is taken into account using the mixing- 
length theory jS ohm- Vitense|[l95^ with the mixing length param- 
eter a = 0. 5. Some exam ples and tests of MAFAGS-ODF models 
are given in'Orupp (2004). Note that due to physical limitations of 
such classical ID LTE models, our analysis is restricted to tumoff 
stars. Only one giant, HD 122563, is included for purposes of com- 
parison with other studies. 

Basic param et ers f or the program stars w ere a dopted 
iGehren et al.l ( |2004 l2006h . and 
These analyses make use of the 



from 



Fuhrmann 



iMashonkina et al 



following methods. The effective temperatures were derived by 
fitting the observed H„ and profiles under LTE assumption. 
The surface gra vities were calculated using parallaxes n measured 
by Hippa rchos (IeSAII 19971) and bolometric corrections were taken 
BC fromlAlonso et al j^ 1995). Ma sses were determined from the 
tracks of iVandenBe rg et alj 1 2000l) by interpolating in the Mbd - 
Jeff diagram. The iron abundances and microturbulent velocities 
were obtained from Fe ii line profile fitting in LTE. For the thick 
disk and halo stars, the abundances of cf-elements are enhanced 
by ~ 0.4 dex. The population identification is based on stellar 
kinematic properties, ages, and abundance ratios [Al/Mg]. All 
stellar parameters are listed in Table[3] 



3.2 Atomic model 

The model of the Co atom is constructed as follows. The num- 
ber of levels and transitions is 246 and 6027 for Co i, and 165 
and 2539 for Co ii. The energy separation of the highest ex- 
cited level from the continuum is 0.4 eV for Co i and ~ 3 eV 



' The differential element abundance in a metal-poor star is given by 
[El/H] = log(^./e). - log(g/£)o 



Table 2. Lines of Co i and Co ii selected for solar and stellar abundance 
calculations. N denotes a number of HFS components for a line, and the 
multiple! is specified in the 2-d column. Ei„„ is the excitation energy of the 
lower level of a transition. refers to the line equivalent width. An asterisk 
in the wavelength entry refers to the lines, which are not used in determina- 
tion of the solar Co abundance due to blending or inaccurate logg/'s. 



A Mult. 

[A] 



N E\av, Lower 
[eV] level 



Upper Wa log gf 
level [mA] 



logCe 



Co I 

3845.470* 34 



3957.930* 

4020.905 

4066.360* 

4110.530 

4121.320 

4792.862* 

4813.480* 

4867.870* 

5212.691 

5280.631 

5301.047 

5331.460 

5352.049 

5369.590* 

5483.340 

5647.234 

6189.000* 



18 

16 

30 

29 

28 

158 

158 

158 

170 

172 

39 

39 

172 

39 

39 

112 

37 



6454.990 174 

6814.950* 54 

7417.380* 89 

7712.661* 126 
Co n 

3501.730 2 



0.92 
0.58 
0.43 
0.92 
1.05 
0.92 
3.25 
3.22 
10 3.17 
9 3.51 
8 3.63 
1.71 
1.78 
3.58 
1.74 
1.71 
2.28 
1.71 
8 3.13 
5 1.96 
8 2.04 
8 2.04 



b'^F, 
a^F. 

' "9/2 
' '-'11/2 

" ^^5/2 

a r|^2 

Z <^ll/2 

" ^^3/2 

" ^^5/2 

" ^^3/2 

" ^^5/2 

Z LI-JI2 

3/2 



,4po 
Z rg^2 

y "^7/2 

5/2 



e^F. 
e^F. 



9/2 
'5/2 



7/2 
9/2 



/'F9/2 



7/2 
o 

5/2 
o 

3/2 
9/2 



y*Fi° 



y '-'7/2 



z*D° 



9/2 



b^P. 



3/2 
3/2 



3/2 

' ^5/2 

7-D° 
^ ^-^5/2 



120 
60 
76 
59 
95 
131 
34 
44 
60 
25 
20 
21 
17 
26 
44 
51 
14 
11 
15 
19 
11 
10 



0.01 2 

-2.07 2 

-2.04 1 

-1.60 1 

-1.08 2 

-0.30 1 

-0.07 3 

0.05 3 

0.23 3 

-0.11 2 

-0.03 2 

-1.94 1 

-1.99 1 

0.06 



2 -: 



-1.59 1 



-1.41 
-1.56 
-2.45 
-0.25 
-1.90 
-2.07 
-1.57 



31.5 
31.7 
31.8 
31.6 
31.6 
31.7 
30.5 
30.5 
30.6 
30.5 
30.4 
31.5 
31.5 
30.5 
31.5 
31.5 
31.4 
31.7 
30.4 
31.7 
31.7 
31.5 



8 2.19 a^Pj 7^D° 86 -1.22 4 -32.2 



" References: (1) iNitz et al 


Jl999l); (2) ICardon et al.l 


iKurucz & PevtremamJ (19961); (4) 


Raassen et al.Hl998l) 



(3) 



for Co II. Energies of levels and wa vele ngths of transitions ar e 
taken from lPickering & Thornj ( Il996h and lPickering et alj (Il998h . 
The oscillator strengt hs are adopted from the Kurucz' database 
jKurucz & Belli 0*9951) . which includes laboratory measurements 
and data calculated using scaled Thomas-Fermi-Dirac radial wave- 
functions. A complete Grotrian diagram for Co is available online. 

Hyperfine splitting of the levels was not included in the SE cal- 
culations, since relative populations of HFS components are ther- 
mal. Moreover, for the uppermost terms above 6.3 eV in Co i the 
fine structure is not maintained. These terms are represented by a 
single level with a weighted mean of statistical weights and ion- 
ization frequencies of their fine structure levels. The transitions 
between two combined levels are also grouped, and the oscillator 
strength of a resulting line is the average of logg/'s weighted ac- 
cording to the appropriate lower-level statistical weights. 

Since no quantum-mechanical calculations of photoionization 
for Co I and Co 11 are available, hydrogenic approximation is 
used. The c ross-sections a re computed from the formula of Kramer 
jMenzel & Pekerislll9"35h corrected for the ion charge (see, e.g., 
lRuttea,2003h and using the effective hydrogen-like main quantum 
number. For cross-sec tions of coUisional ionization, we use the for- 
mula of lSeatonl ( fT963) . The rates of allowed and forbidden bound- 
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Table 3. Stellar parameters and their estimated errors for the selected sample. In the Ref column, the sources of adopted stellar parameters are given. 
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HIP 


Instr. 
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[Fe/H] 


[Mg/Fe] 




Population 


n J-' (3 

Refj^ 


Nco 


[Co/F 


e] 








[K] 




[km/s] 














NUTE 


LTE 


HD 19445 


14594 


FOCES 


5985 ± 80 


4.39 ±0.05 


1.5 


-1.96 


0.38 


22.7 


Halo 


1 


1 


0.64 


0.08 


HD 25329 


18915 


FOCES 


4800 ± 80 


4.66 ±0.08 


0.6 


-1.84 


0.42 


50.1 


Thick Disk? 


4 


1 


-0.02 


-0.72 


HD 29907 


21609 


UVES 


5573 ± 100 


4.59 ±0.09 


0.9 


-1.60 


0.43 


17.3 


Halo? 


2 


2 


0.66 ± 0.06 


0.19 ±0.05 


HD 34328 


24316 


UVES 


5955 ± 70 


4.47 ± 0.07 


1.3 


-1.66 


0.42 


14.4 


Halo 


2 


3 


0.67 ± 0.01 


0.11 ±0.16 


HD 61421 


37279 


UVES 


6510 ± 100 


3.96 ±0.05 


1.8 


-0.03 


0.0 


324.9 


Thin Disk 


5 


8 


-0.09 ± 0.07 


-0.15 ±0.1 


HD 84937 


48152 


UVES 


6346 ± 100 


4.00 ±0.08 


1.8 


-2.16 


0.32 


11.7 


Halo 


1 


3 


0.66 ± 0.03 


0.18 ±0.01 


HD 102200 


57360 


UVES 


6120 ±90 


4.17 ±0.09 


1.4 


-1.28 


0.34 


10.5 


Halo 


2 


3 


0.51 ±0.06 


0.03 ±0.11 


HD 103095 


57939 


FOCES 


5110± 100 


4.69±0.10 


1.0 


-1.35 


0.26 


140.0 


Halo 


1 


1 


-0.01 


-0.11 


BD-4°3208 


59109 


UVES 


6310 ±60 


3.98 ±0.21 


1.5 


-2.23 


0.34 


3.7 


Halo 


2 


2 


0.71 ±0.08 


0.13 ±0.01 


HD 122196 


68464 


UVES 


5957 ± 80 


3.84±0.11 


1.7 


-1.78 


0.24 


7.4 


Halo 


2 


2 


0.52 ± 0.02 - 


-0.12 ±0.06 


HD 122563 


68594 


UVES 


4600 ± 120 


1.50 ±0.20 


1.9 


-2.51 


0.45 


5.2 


Halo 


3 


3 


0.57 ± 0.02 - 


-0.07 ± 0.06 


HD 134169 


74079 


FOCES 


5930 ± 200 


3.98 ±0.1 


1.8 


-0.86 


0.53 


15.1 


Thick disk 


1 


3 


0.45 ± 0.02 


0.18 ±0.14 


HD 148816 


80837 


FOCES 


5880 ± 200 


4.07 ± 0.07 


1.2 


-0.78 


0.36 


27 


Thick disk 


2 


6 


0.38 ± 0.03 


0.19 ±0.06 


HD 184448 


96077 


FOCES 


5765 ± 200 


4. 16 ±0.07 


1.2 


-0.43 


0.47 


30.4 


Thick disk 


2 


11 


0.34 ± 0.04 


0.27 ± 0.03 


HD 140283 


76976 


UVES 


5773 ± 60 


3.66 ±0.05 


1.5 


-2.38 


0.43 


18.0 


Halo 


2 


3 


0.85 ± 0.02 


0.14 ±0.01 


G 20-8 


86443 


FOCES 


6115 ±80 


4.20 ±0.20 


1.5 


-2.19 


0.45 


5.1 


Halo 


1 


2 


0.76 ± 0.06 


0.28 ±0.15 


G 64-12 


66673 


UVES 


6407 ± 80 


4.20 ± 0.20 


2.3 


-3.12 


0.33 


1.3 


Halo 


1 


1 


0.8 ±0.1 


0.19±0.1 


HD 200580 


103987 


FOCES 


5940 ± 80 


3.96 ±0.06 


1.4 


-0.82 


0.46 


13.8 


Thick Disk 


2 


2 


0.38 ± 0.04 


0.12 ±0.1 



" References: (n iGehren et"ai] b006l) ; (2) lGehren et all j2004h ; (3') lMashonkina et al.l i2008l) ; (4) lFuhrmanr] i2004 ; fS^ lKom et al.l i2003h 



bound transitions du e to collisions wi t h elec trons are ca l culate d us- 
ing the formulae of Ivan RegemorteJ ( Il962h and I Alleril ( Il973h . re- 
spectively. 

The rates of bound-bound and bound-free transitions due to in- 
elastic collisio ns with H i atom s are computed from the formula of 
lDrawinl i [l969l) in the version o f ISteenbock & Holwegeil h984h and 
are multiplied by a scaling factor 5h = 0.05. This value gave the 
smallest fitted-abundanc e spread in our previous work on SE of Mn 
in stellar atmospheres ( Bergemann & Gehren 2007, 2008). Most 
abundance analyses favour simila r small values of 5h (M ashonkina 
1994 lAbia & Mashonkinall2004 ISundgvist et ai]|2008l : IShi et al.l 
20081) . To the best of our knowledge, there are onl y a few indi- 
cations in the literatur e that 5h is larger than unity jOratton et al.l 
ll999 ^:'Kometal.'2003). Unfortunately, these cases refer to the most 
prominent member of the Fe-group, iron. 

In subsequent sections, we study the influence of photoioniza- 
tion cross-sections and inelastic collisions with hydrogen on popu- 
lations of atomic levels in Co. However, it should be kept in mind 
that values for cross-sections have only order-of-magnitude accu- 
racy. 



4 RESULTS 

4.1 Statistical equilibrium of Co 



To discuss the excitation-ionization equilibrium of Co in stellar 
atmospheres, we use the notation of the departure coefficient hj, 
which is defined as bj = N^™/N™ with A', the number density of 
atoms in the level The ^),-factors as a function of optical depth at 
500 nm are presented in Fig.|2]for the solar model atmosphere and 
several stellar models from the grid. A few selected levels, typical 
for their depth dependence, are indicated by heavier black traces; 



the remaining bulk of the Co i and Co n levels are indicated by grey 
lines. 



4.1.1 The Sun 

Fig.|2ji shows the departure coefficients, calculated with our refer- 
ence model of Co for the Sun. The overall underpopulation of all 
Co I levels is due to the overionization mainly from the levels of 
2 - 3 eV excitation, e.g. z~G° with threshold wavelength 3151A. 
The Co I ground state is not subject to strong overionization. But 
its separation from the first metastable level is only 0.4 eV, which 
favours a very efficient coUisional interaction. Strong collisional 
and radiative coupling also persists between other Co i levels, be- 
cause the total number of transitions in Co i is very large. Thus, 
the distribution of ii-factors up to logrjooo ~ -2 has a very reg- 
ular structure, with levels of ever decreasing energy gap from the 
continuum showing smaller deviations of from unity. 

In the higher layers, spontaneous transitions in strong lines 
of Co I slightly modify this simple pattern. At logTsooo < -1.5, the 
medium becomes transparent to the radiation in transitions between 
the low metastable levels and the levels of odd doublet and quartet 
term systems with ~ 4 eV excitation energy, that leads to a deple- 
tion of the upper levels. One such transition, a^F^^j — > z'Ggp, is of 
particular interest, because the corresponding spectral line 4121 A 
is the only observable line of Co i in our spectra of very metal-poor 
stars. 

A large number of strong transitions with UV wavelengths in 
Co I connect the low levels with - 3 eV and the levels with 5 - 6 eV 
excitation energies Hence, at depths with -1 < logTsooo < 0.2 
there is also overexcitation of the upper levels. This maintains 
nearly constant Zj, of upper levels at - 1 < log Tjooo < 0, in spite of 
increasing overionization (e.g. level e in Fig. |2^). This process 




Figure 2. Departure coefficients fo, of selected Co i and Co ii levels for the solar model and for several stellar models from the grid. 



barely affects populations of the lower levels due to the small Boltz- 
mann factor, and dominance of photoionization from these levels. 

The majority of the uppermost levels with £, > 6.5 eV 
are underpopulated relative to the continuum. They interact much 
stronger radiatively with the low levels and with each other by 
means of collisions than with the Co ii ground state, a^F^. Only 
for a few highest levels, with threshold ionization energies similar 
to the mean kinetic energy of the electrons in the solar atmosphere, 
does the collisional coupling to the a 'F^ level prevail. Thus, they 
are in thermal equilibrium with a ^F^. 

The ground state of Co ii is also affected by NLTE, although 
departures in the line formation region are not large. Low-excitation 
levels of Co n are in thermal equilibrium with its ground state. 
But, most of the higher Co n levels are overpopulated due to non- 
equilibrium excitation processes. This NLTE mechanism is typical 
for majority ions, i.e. ions of the dominant ionization stage of an 
element. 

The tests performed with a smaller model constructed with 
246 Co I levels and closed by the Co n ground state indicate poor 
performance of such a reduced atomic model. This changes ioniza- 
tion and excitation equilibria in both ions, which is manifested in 
spurious overpopulation of the Co ii ground state and stronger ther- 
malization of Co i levels. Thus, reduction of the atomic model, as 
is usually done for 3D NLTE calculations, must be performed with 
extreme caution. 

The uncertain parameters in our models are cross-sections 
for photoionization and collisions with neutral hydrogen atoms. 
Fig. I2J3 demonstrates the effect of stronger photoionization on the 
atomic level populations; here, cross-sections are multipled by a 
scaling factor Sp = 300. The main difference with the reference 
model (Fig. [2^) is that depopulation of all levels is amplified at 
-1 < logTsooo < 0.4. The influence of collisions with H i is fully 
understood from Fig. [2j;. When the standard Drawins formula is 
used with the scaling factor 5h = 1, all levels at line formation 
depths -2 > logTsooo > are more strongly coupled to the con- 
tinuum, but this scaling factor is too small for LTE to be reached 
completely. 



4.1.2 Metal-poor stars 

The character of processes leading to NLTE effects in Co does not 
change for stellar atmospheres with parameters of cool turnoff stars, 
which form the basis of our abundance analysis. 

Fig. |2ji demonstrates the effect of increasing temperature on 
the departure coefficients of selected Co i and Co 11 levels. At 
Teff = 6400 K, the ionization equilibrium Co i/Co 11 is somewhat 
different from that of the solar model. Nearly 98% of Co atoms are 
in a singly ionized stage, hence the ground state of Co 11 is not sen- 
sitive to overionization from the Co i levels. The level a^F^ keeps 
its thermodynamic equilibrium value. However, departure coeffi- 
cients in Co I show larger deviations from unity. The reason is that 
a stellar flux maximum is shifted to shorter wavelengths, thus ove- 
rionization affects even the lowest metastable levels with £, ~ 1 eV 
and /liiu ~ 1800 A. In contrast, in the solar model overionization is 
important only for the levels with £, ~ 3 eV. 

The influence of photoionization is more pronounced with de- 
creasing metallicity (Fig. [2^). The dramatic depopulation of all lev- 
els occurs at the depths of line formation, -1 < logTjooo < 0. 
The overionization from the lowest levels is balanced by increased 
net recombination to the higher levels that results in a small over- 
population of the latter at logTsooo ~ 0. In the higher layers, 
logTsooo ~ -1 even an overpopulation of the Co 11 ground state 
develops. Also, we note an increasing importance of line pump- 
ing in Co II. For example, the level z^DJj has b, » I already at 

log T5000 0.3. This is characteristic of all Co 11 levels with £, > 5 

eV. Hence, we expect substantial NLTE effects in the lines of both 
Co I and Co 11 in metal-poor stars. 

A combination of low effective temperature, low gravity, and 
low metal content leads to a rather disorderly behavior of depar- 
ture coefficients for different levels (Fig.|2f). Such stellar parame- 
ters are representative of giants where the collisional interaction is 
much weaker compared to dwarfs. Each level of Co i is no longer 
coupled to the bulk of levels with similar energies, but displays a 
distinct behaviour that is determined by the competition of radiative 
b-f and b-b processes. As a result, NLTE line formation for giants 
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Figure 3. The lines of Co i (4121, 6082 A) and Co ii (3501 A) in the solar 
flux spectrum (filled circles). Synthetic NLTE and LTE profiles with and 
without HFS are labeled correspondingly. Wavelength positions and relative 
line strengths of the HFS components are indicated. Note that the relative 
HFS component intensities are not on the same scale in figures. The long- 
dashed trace on the bottom figure indicates the NLTE profile computed with 
A = 49 mK for the lower level a ^P, (our reference value is A = 40 mK). 



and dwarfs can be different for transitions involving different lev- 
els. This is evident by comparing the runs of bi for the levels of 
e^'G, z^G°/2, and a'F terms in Fig.|2t and|2f. 



4.2 NLTE and HFS effects on the fine formation 

The computed profiles for selected fines of Co i and Co ii are com- 
pared with the solar spectrum in Fig. [3] The synthetic profiles were 
calculated under LTE and NLTE conditions, with and without HFS. 



Table 4. NLTE abundance corrections for lines of Co i and Co n calculated 
with selected models of the grid. Hyphens refer to the fines with computed 
NLTE equivalent widths below 1 mA. Colhsions with neutral hydrogen are 
included with the reference scaling factor Sh = 0.05. Note that Anlte for 
the Co II line at 3501 A are also given for the cases when its equivalent 
width is below 1 mA. 



reff/logg/[Fe/H] 
wavelength, A 



Anlte 
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All lines selected for the abundance analysis in the Sun and 
metal-poor stars are relatively weak (Table|2j. Hence, the Co i lines 
computed under NLTE behave similarly: they are uniformly weak- 
ened compared to the LTE profiles. The main effect is the shift of 
the optical depth scale due to fo, < 1 at line formation depths. As a 
result, the NLTE abundance correction^ Anlte are positive; for the 
majority of Co i lines computed with the solar model atmosphere, 
Anlte — hO.15 dex. 

NLTE abundance corrections for different stellar parameters 
are given in Table |4] These values are computed only for five Co i 
lines and one Co ii line, which can be detected in our spectra of 
stars with [Fe/H]< -0.5. Although it is not our goal to investigate 
giant stars, we have performed test calculations for the cool giant 
model {T^s = 4800 K, logg = 1.8, [Fe/H] = -3.3). This exception 
is made because many analyses of Co abundances refer to bright 
cool giants with low [Fe/H]. NLTE corrections from Table |4] should 
only be used to get a crude estimate of deviations from LTE in 
metal-poor stars due to the dependence of NLTE correction on the 
Co abundance. 

The main stellar parameter that determines the sign and mag- 
nitude of the NLTE abundance corrections is metallicity. Low metal 
abundances lead to decreasing number density of free electrons and 



The difference in abundances required to fit NLTE and LTE profiles is re- 



ferred to as the NLTE (abundance) correction Anlte = log e 
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Figure 4. NLTE abundance corrections Anlte calculated for the lines of 
Co I (4121 A, top) and Co n (3501 A, bottom). Calculations are performed 
for sixteen models with T^s = 5000, 5500,6000 K, logg = 3,4, [Fe/H] = 
0,-1,-2,-3. 

reduced line absorption in the UV, so that the effect of overioniza- 
tion on the opacity of Co i hnes monotonously increases. On av- 
erage, Anlte increases with decreasing [Fe/H] from ~ 0. 1 dex for 
[Fe/H] = to ~ 0.7 dex for [Fe/H] = -3 (Fig.|l]top). It is clear that 
in addition to the effect of metallicity, high temperatures control 
overionization at higher [Fe/H], whereas at low [Fe/H] the effect of 
low log g is more important. 

The response of NLTE corrections for Co i to the effective 
temperature is not uniform. At moderate gravities (log g = 4), the 
NLTE corrections are maximal in the warm model T^g = 5500 
with the smallest metallicity (see Table |4l(. Lower NLTE correc- 
tions at even higher temperatures are due to the increased rates of 
collisions. At low temperatures, NLTE abundance corrections are 
also smaller because the stellar flux maximum is shifted to longer 
wavelengths, away from ionization thresholds of the important Co i 
levels. 

NLTE corrections increase with decreasing surface gravity. 
However, for different metallicities and supersolar temperatures, 
reduction of logg from 4.6 to 3.4 changes the NLTE abundance 
corrections for all lines only by ~ 0.03 - 0.06 dex. In contrast, 
Anlte strongly depends on stellar gravity in the coo! model with 
Teff < 5000. The extreme values of Anlte are found for the cool 
giant model with a subsolar temperature T^g = 4800 K, [Fe/H] 
= -3.3, and a very low gravity logg = 1.8. This is an expected 
result, because all collisional interactions become very weak due to 
the small metallicity and gravity, and the role of radiative processes 
increases in spite of the reduced temperature. 

The NLTE corrections for the Co ii line show a different be- 
haviour, which is related to the deviation of line source function 
S' from the Planck function By{T^). As an example, we study the 
formation of the Co ii line at 3501 A. It seems to be relatively un- 
blended in the solar spectrum and can be used to test the ionization 
equilibrium of Co in the metal-poor stars. In the solar model, this 
line is sufficiently strong; its core is formed in the upper layers 

log T5000 2, where spontaneous transitions depopulate the upper 

level, thus bi > bj and 5' < By{T^). Therefore, line core inten- 
sities are slightly larger in NLTE (Fig. [3] bottom), and Anlte < 
(Fig.|4]bottom). At [Fe/H] < -1, NLTE abundance corrections start 
to increase smoothly. Due to less UV opacity, the radiation field 



is amplified producing overpopulation of the higher levels via the 
photon pumping mechanism, b, < bj. The line source function be- 
comes superthermal, S ' > By(T^), and the line is weakened com- 
pared to the LTE case. The same mechanism is responsible for the 
positive NLTE abundance corrections in the models with [Fe/H] 
< -2. At very low metallicities (Fig.|4]bottom), NLTE corrections 
are as large as +0.2 . . . + 0.3 dex, depending on the temperature. 
Consequently, the analysis of Co 11 lines in very metal-poor stars 
can lead to substantial errors in abundances if NLTE effects are 
neglected. 

The errors in abundances introduced by neglecting HFS de- 
pend on the line strength. For the solar model, the abundances de- 
rived for saturated lines 4I2I A and 3501 A without HFS are over- 
estimated by 0.4 dex and 0.9 dex, respectively. The reason is that 
HFS effectively de-saturates strong lines leading to a compound 
profile, where the strength of each component is more linearly pro- 
portional to the element abundance. Note also that saturated lines 
are also sensitive to the A and B factors used in calculations of 
HFS. For example, the 3501 A HFS line computed with an approx- 
imate estimate A = 49 mK for the lower level a ^P, (taken from 
[Pickering et alj[l99^ cannot fit the observed profile at all (Fig. [3] 
bottom). In contrast, our new reference value A = 40 mK provides 
an excellent fit to the observed profile. The weak lines, like 6082 
A, are less affected: the profile computed with all HFS components 
gives ~ 0.1 dex lower abundance. At low metallicity, [Fe/H] = -2, 
the abundances are almost insensitive to inclusion of HFS in spec- 
trum synthesis; the errors for all investigated lines are not larger 
than 0.05 dex. 

4.3 Solar abundance of Co 

The solar spectrum is calculated using the MAFAGS-ODF model 
atmosphere with T^tf = 5780 K, logg = 4.44, [Fe/H] = 0, and a 
constant microturbulence velocity = 0.9 km s ' for all lines. The 
Co line profiles are broadened by a rotation velocity Vmt = 1.8 km 
s"' , and by a macroturbulence velocity V^.^c = 2.5 ... 4 km s"' . The 
observed flux spectr um is taken from the Kitt Peak Solar Flux Atlas 
( iKurucz etaljl984h . 

The lines of Co i and Co 11 selected for abundance calculations 
are listed in Table [2] Measured oscillator strengths are available 
only for 17 out of 20 lines. Two refere nce s ources of experim ental 
logg/ values are lCardon et al.l ( Il982h and iNitz et al] ( Il999 ). For 
three yello w lines of multiplet 1 5 8 we had to adopt the theoretical 
values of Kuruc z & PevtremannI ( Il99d ). However, the accuracy of 
calculated data is unknown, hence the lines of multiplet 158 are 
only used in the differential abundance analysis of metal-poor stars. 
Our final estimate of the solar Co abundance is based solely on the 
experimental log gf values. Three lines in the near-UV, 3845, 3957, 
and 4066 A, are also neglected in solar calculations due to severe 
blending. 

Figure |5] shows NLTE and LTE abundances for 16 lines of 
Co I as a function of their oscillator strengths. A weak trend of in- 
dividual abundances with gf values, as well as similar line-to-line 
scatter under NLTE and LTE, points at th e errors in logg/' s. For 
the weak lines, the errors are ~ 15-25% l iCardon etaljl982h . Two 
weak lines at 6189 and 6814 A show NLTE and LTE abundances 
deviating from the mean by > 0.1 dex. The lines at 7417 and 7712 
A contain a well-resolved blend in their wings, which is, however, 
not present in our linelist. The line at 5369 A gives very large abun- 
dance, 5.15 dex. As its g/-value seems to be reliable, it is likely that 
the line contains an unidentified blend. We ignore these five lines in 
the solar analysis, and our sample is further reduced to 1 1 lines. No 
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significant correlation of individual NLTE or LTE abundances is 
found with the excitation potential of the lower level or equivalent 
width of the lines. 

The average NLTE and LTE abundances of Co in the solar at- 
mosphere are log = 4.95±0.04 dex and log ej;™ =4.81 ±0.05 
dex, respectively. The standard deviations are quot ed as errors. 
When calculated only with the log gf data from ICardon et al.l 
('l982), log = 4.96+0.06 dex. This value is nearly equal to the 
NLT E abundance calc ulated using 5 lines with oscillator strengths 
from lNitzetaDjl999l) . loge^™ = 4.98 ± 0.03. 

It is interesting that there is a large difference of 0. 14 dex be- 
tween the solar NLTE and LTE abundances determined from the 
Co I lines with Sn = 0.05. Up to now, such a strong NLTE abun- 
dance effect has been demon strated only for Sc i, with Anlte ~ 0. 18 
dex derived using Sh = 0.1 jZhang et alj|2008[) . The average solar 
NLTE corrections for the other Fe-peak elements, Fe jGehren et al.l 
1200 ih and Mn (Bergemann & Gehren 2007), are smaller. This is, 
however, not surprising because the magnitude of NLTE effects 
strongly depends on the model atom used in calculations, in partic- 
ular, on the collision efficiency with hydrogen atoms 5h (see Sect. 
|3.2| l and the detailed atomic properties, like separation of energy 
levels, number and strength of transitions, and ionization potential. 
Weaker NLTE effects for Mn i (Anlte ~ 0.05 dex) compared to 
Co I are due to the differences in the atomic models: a) the num- 
ber of low-excitation levels, which stipulate overionization; b) col- 
lisional co upling to the c ontinuum (the latter issue was also em- 
phasized in lOehren et al.ll2001) , which amplifies recombination; c) 
the general number of levels. The model of the Fe i atom used by 
lOehren et alj ( 1200 1[) is different from our Co i model (at least) in 
the following aspects. The authors adopted 5h = 5 as their final 
value and forced thermalization of the Fe i levels with excitation 
energy above 7.3 eV. Assuming no explic it thermalization, as we do 
in the current work, iGehren et alj ( 1200 ih obtain very strong NLTE 
effects on Fe i lines even for the model with 5h = 5, Anlte ~ 0. 12 
dex. Using this large scaling factor for H i collisions, we derive the 
average NLTE abundance correction for Co i lines Anlte ~ 0.06 
dex. Thus, we see no reason to force thermalization of the upper 
Co I levels. In contrast, we consider collisions with H i as the main 
source of log uncertainty. We performed a series of NLTE 
abundance calculations with < 5h < 5. The NLTE result for 
Sh = 0, loge^™ = 4.96 ± 0.04, is almost identical to the case of 
Su = 0.05. Also, with increasing scaling factor from 0.5 to 5, the 
abundance steadily decreases from 4.91 to 4.87 dex with the stan- 
dard deviation remaining constant, cr » 0.04 dex. Thus, no con- 
clusion can be drawn concerning the optimum value of the scaling 
factor to hydrogen collisions. 

The majority of Co i lines are weak, and they are not sen- 
sitive to the van der Waals damping, neither are they affected by 
the microturbulence. We found only one line at 4121 A, which 
gives different abundances with the variation of or logCe. For 
A^, = +0.2 km s"', Alog£ = +0.1 dex; AlogCe = +0.5 corre- 
spond s to Alogg = +0. 12 dex. With logCe = -31.63 calculated 
using lAnstee & 0'Mm3 (fl995) theory and ^, = 0.9 km s"', this 
line gives a NLTE abundance of 4.95 dex, which is consistent with 
abundances derived from the other Co i lines. Moreover, the stan- 
dard deviation of the average Co NLTE abundance, cr = 0.04 dex, 
is already smaller than the errors given for the oscillator strengths. 
Hence, we have not tried to adjust damping constants and micro- 
turbulence any further. 

We checked the ionization equilibrium of Co for the Sun, us- 
ing the relatively unblended line of ionized cobalt Co ii at 3501 A 
(Fig- m bottom). The line is saturated and consists of 8 HFS com- 
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Figure 5. Abundances of tlie 17 solar Co i lines. NLTE and LTE abundances 
are marked with filled and open circles, respectivety. 



ponents, which were calculated with our new hyperfine splitting 
constants. To produce a better fit to the observed profile, the A fac- 
tor for the upper level was decreased by 1 mK; this is within the 
errors of the me asured HFS value. U sing the theoretical oscilla- 
tor strength from lRaassen et af] ( Il998t) . log gf = -1.22, we derive 
loge^™ = 4. 86 and log = 4.88 dex. With the experimental 
g/- value from ISalihet all (1 19851) 



1), logg/ 



-1.18, the NLTE and 
LTE abundances are loge^™ = 4.82 and log£^™ = 4.84 dex, 
respectively. We do not assign any error to these values for the rea- 
sons described in the following paragraph. 

The ~ 0. 1 dex discrepancy between the lines of Co i and Co ii 
in NLTE is most likely due to the erroneous abundance from the 
Co II line. First, the uncertainty due to the continuum placement 
and blending in the near-UV window of the solar spectrum is large. 
At present, there is no information about other strong lines overlap- 
ping with the 3501 A line of Co ii. However, if atomic physics ex- 
periments or theoretical calculations reveal the presence of a strong 
blend, the Co abundance determined from the UV Co ii line will be 
re duced. The accuracy o f the calculated oscillator strength provided 
bvlRaassenetal.l ( ll998[) is not known, whereas the accuracy of the 
ISalih et all ( Il985h g/-value for the 3501 A line is not better than 
50%, which translates to abundance uncertainty of ~ 0.2 dex. If the 
problem is concealed in hydrogen collisions, a very large scaling 
factor 5h > 5 is needed to bring two ionization stages in agreement. 
This issue will be further investigated in the next section devoted to 
the metal-poor stars. 

The abundance of Co in C i meteorites is 4.89 + 0.01 
dex dUodders et al. I l2009l) . It conforms to our NLTE abundance 
loge^™ = 4.95 ± 0.04 dex derived from Co i lines within the 
combined errors of both values. Also, the LTE abundance from 
the Co II line (4.88 dex) agrees to meteoritic, if the theoretical 
g/-value from Raassen et al. (1998) is used. Previously, the pho- 
tospheric LTE ab undance from the Co i lines was determined by 
ICardonetal]fl982i) . logEco.o = 4.92±0.08 dex. This value is ~ 0.1 
dex higher than our LTE abundance, which is due to several rea- 
sons. First, the authors used the semi-empirical model atmosphere 
of iHolw eger ( 1 967), which is the essentially the same as the model 
of iHolweger & Mtilleil ( Il974l) . The latter model is known to give 
higher abundances compared to the theoretical model atmospheres. 
Second, equivalent widths were used to measure the abundances. 
This method always gives only an upper limit on the abundance. 
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Figure 6. [Co/Fe] vs. [Fe/H] in metal-poor stars. The Co abundances are 
derived from Co i lines under NLTE (filled circles) and LTE (open circles). 
Iron abundances are determi ned from Fe n lines in LTE. Th e sources of 
other LTE [Co/Fe] data are: Johnson"2002' (open triangles). iReddv et alj 
I2003ll2006l (crosses). [McWilliam et al...l995. (stars). 

because the contribution to a line strength from blending lines is ne- 
glected. In all other respects (oscillat or strengths, microtu rbulence, 
line set), our study and the analysis of ICardon et alj l [T982l) are sim- 
ilar. 



4.4 Abundances in metal-poor stars 

Our sample of 17 objects from iBergemannI ( |2008|) was expanded 
by one subdwarf G 64-12 with [Fe/H] = -3.12. The UVES spec- 
trum for G 64-12 was taken from ESO/ST-ECF Science Archive 
Facility (PID 67.D-0554(A)). Other stars were observed with the 
UVES echelle spectrograph (ESO VLT UT2, Chile) in 2001, and/or 
with the FOCES echelle spectrograph (CAHA observatory, Calar 
Alto) during 1999 and 2000. The data fo r HP 61421 and HP 8 4937 
were taken from the UVESPOP survey (Bagnulo et al ] |2003l) . The 
UVES data cover a spectral range froin 3 300 to 6 700 A, whereas 
FOCES spectra are cut at 4000 A. A de tailed description o f the 
observational material can be found in iBergemMii] ^20081) and 
iBergemann & GehrerJ |2003) . 

The abundance ratios for Co with respect to Fe in stars of our 
sample are given in Table[3] The Co abundances of thick disk stars 
with -1 < [Fe/H] < -0 are based on 6 to 10 lines, which have W.i > 
20 mA in the solar spectrum. In the spectra of the halo stars, only a 
few lines of Co i {AA 3845, 4020, 41 10, 4121) are detected. Pue to 
a systematic abundance offset of ~ 0.2 dex between the line at 41 10 
A and other Co i lines, this near-UV line was excluded from stellar 
analysis. The discrepancy might be related to the poor modelling of 
Ha inner wings, where the 41 10 A line is located, or to a blend. In 
fact, no Co abundance studies in metal-poor stars include the line at 
41 10 A. For G 64-12, the abundance is based on the single Co i line 
at 3845 A, and the error bar ±0. 1 dex reflects the uncertainty of the 
line fitting due to the low S/N (S/N ~ 120) in the near-UV spectral 
window. We could not determine the solar log(g/£)Q value for the 
3845 A line, hence all stellar abundances based on the line were 
calcul ated using the experimental logg/ = 0.01 from lCardon et al.l 
( Il982l) . 

The [Co/Fe] - [Fe/H] relation, along with the results of other 
authors, is demonstrated in Fig. [6] Note that the solar NLTE and 



LTE [Co/Fe] values are set to by definition. The NLTE abun- 
dances of Co in the metal-poor stars are higher than the LTE abun- 
dances by 0.4 - 0.7 dex. Thus, instead of the flat [Co/Fe] trend 
in metal-poor stars, we find that the [Co/Fe] ratio increases with 
decreasing metallicity down to [Fe/H] ~ -2 and is roughly con- 
stant thereafter. Two stars with very low NLTE and LTE [Co/Fe] 
ratios a re HP 25 329 and HP 103095. These stars were also anal- 
ysed by Peterson (1981), who derived [Co/Fe] = +0.48 and [Co/Fe] 
= -1-0.15 for HP 25329 and HP 103095, respectively. The discrep- 
ancy with our values stems from dilTerent model parameters and 
abundance analysis techniques. In particular, Peterson used log gf 
values calibrated on solar lines to calculate the abundances from 
line equivalent widths, and Unsold approximation for log Ce values 
was assumed. 

Our LTE abundances are generally consistent with LTE results 
of the other authors considered here, although in the metal-poor 
domain Co abundances show a large scatter within a study and be- 
tween studies (Fig.|6j. The general trend is hard to define, although 
the majority of studies indicate that Co is increasing with respect 
to iron towards very low metallicities. For the thick disk stars our 
[Co/Fe] trend with [Fe/H] is flat and slightly supersolar, in agree- 
ment with Reddy et al . (2006). The thin disk star, Procyon, supports 
mild un derabundance of Co in the thin disk as follows from the re- 
sults of iReddv et all ( |2003|) . The average Co abundances for stars 
[Fe/H] < -1 are consistent with ljohnsoni ( |2002|) as to slope, abso- 
lute values, and dispersion. Note that Johnson analyzed only giants. 

In the attempt to check the ionization equilibrium of Co, we 
calculated differential abundances from the 3501 A line of Co n. 
This line can still be discerned from the continuum only in the 
spectra of 5 stars. Thus, these abundances must be treated, at best, 
as upper limits because of poor S/N below 4000 A. The NLTE 
abundance corrections for Co n are not so small andpositive. For 
HP 84937 and HP 140283, A[Co/Fe]„(NLTE - LTEf^ -1-0.14 and 
+0.11 dex, respectively. Whereas, for HP 102200, HP 34328, and 
HP 122563, the NLTE corrections are +0.04 . . . + 0.07 dex. A sig- 
nificant discrepancy still remains between abundances determined 
from the lines of Co i and Co ii. 

Consistent ionization equilibria under NLTE require either 
significantly lower abundances from the neutral lines and/or much 
higher abundances from the ionic line. The former can be achieved 
with increasing the scaling factor to the H i collision cross-sections, 
5h. Test calculations were performed for two stars. The NLTE re- 
sults for HP 34328 computed with Sh = 1 are not very different 
from the case of 5h = 0.05 (our reference value): the abundance 
from the Co i lines decreased by 0.15 dex. It is only with 5h = 10 
that the discrepancy is solved: [Co/Fe]J^™ = 0.26, which is in 
agreement with [Co/Fe] = 0.21. The abundances derived from 
Co n lines are not sensitive to a variation of 5h. However, for HP 
84937, assumption of 5h = 10 leads to [Co/Fe]i = 0.27, which is 
still « 0.2 dex higher than the NLTE abundance determined from 
Co n lines. Whether such a large scaling factor is justified is ques- 
tionable. 

Another option is that our model atmospheres are not satis- 
factory, and granular inhomogeneities will have a large effect on 
abundances in addition to NLTE. Let us consider Fe, for which cal- 
culations of 1.5P NLTE line formation with 3P convective model 
atmospheres have been performed for the Sun and the metal-poor 
subgiant HP 140283. Using the Fe i fines .Shchukina et al.l ( |2003) 



' Co abundance ratios determined from the Co i and Co n lines are denoted 
as [Co/Fe]i and [Co/Fe]n. respectively 
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found equal 3D and ID NLTE abundances for both stars that were 
interpreted as due to stronger NLTE effects in 3D than in ID mod- 
els, which exactly cancel large 3D corrections. However, large posi- 
tive NLTE abundance corrections of — 1-0.4 dex were demonstrated 
at low metallicity for weak Fe n lines in 3D. Recalling that Co i and 
Fe I have a very similar complex atomic structure with ionization 
energies differing by ~ 0.04 eV, one may expect that 3D will have 
the same (minor) effect on NLTE-based abundances of Co deter- 
mined from Co I lines. Whereas, large positive 3D corrections may 
be necessary for Co n lines. In this case, our metallicities based on 
LTE abundances from Fe ii lines may not necessarily be wrong, 
because the magnitude of the NLTE effect depends strongly on 
the choice of the scaling factor Sw to the hydrogen collision cross- 
sections. This parameter must be separately calibrated on the solar 
and stellar spectra, and for Fe it may take the values different from 
that used for Co. 

We have also investigated the effect of a temperature enhance- 
ment on the abundance analysis by repeating the calculations for 
HD 102200 with Areff = ±100 K with respect to our reference 
value TefF = 6120 K. For both ions the effect turned out to be neg- 
ligible. With = 6220 K, abundances for Co i lines increased 
respectively by 0.02 and 0.06 dex, and the abundances from the 
Co n line decreased by ~ 0.02 dex. Abundance corrections of the 
same magnitude but with the reverse sign were derived for lower 
T^g. Thus, we cannot attribute the ionization equilibrium problem 
to systematic errors in T^ff. 



see that all chemical evolution models, which utilize metallicity- 
dependent yields, are inadequate to describe the [Co/Fe] trend sug- 
gesting that the problem is in the stellar yields. This is not un- 
expected, gi ven the sensitivity o f super nova yields to the details 
of explosion. IWooslev & WeaveJ fl99l) emphasized that their SN 
II iron yields could be in error by a factor of 2, ascribing even 
larger uncertainties to the elements produced in a-rich freeze-out 
due to the great sensitivity of those to the mass cut placement. For 
deeper mass cuts, more Co is produced. iNakamura et"ai] ( 1 19991) 
foun d that an increasing neutron exc ess increases the Co produc- 
tion. ' Umeda & Nomotj ( l2003l I2OO5I) showed that if a large-scale 
uniform mixing occurs in the inner layers of very massive stars, 
M > 25M0, then Co is produced more efficiently and the effect 
is e ven enhanced fo r large r explosion energies. In connection with 
this. lNagataki et al] ( Il997h demonstrated that the axisymmetric ex- 
plosion in SN II models also boosts an tt-rich freeze-out, and hence 
an increased production of Co. Interestingly, the results of massive 
star nu cleosynthesis calculations performed by iHeger & WooslevI 
1 I2OO8I) suggest that there is no need for hypemova to explain ele- 
ment abundance patterns in extremely metal-poor stars. The vari- 
ety of mechanisms, which could affect the Fe-peak production, is 
overwhelming (Thielemann et al. 2007). Thus it is not hard to be- 
lieve that the disagreement between observations and predictions of 
chemical evolution models is, at least in part, due to the deficiencies 
of the latter. 



5 COMPARISON WITH CHEMICAL EVOLUTION 
MODELS 

There are several models in the literature that describe the evo- 
lution of Co in the Galaxy. The difference between theoretical 
trends is strikingly large: both the shape of the [Co/Fe] trend and 
the magnitude of Co depletion or overproduction relative to Fe 
in metal-poor stars vary by an order of magnitude. In fact, nei- 
ther of the GCE models are able to describe our NLTE [ Co/Fe] 
trend in metal-poor stars. The lGoswami & PrantzosI ( I2OO0I) model 
with metallicity-independent SNe yields gives only a qualitatively 
similar behaviour, characterized by increase of [Co/Fe] ratios to- 
wards low metallicities. However, this model gives a very good 
fit to the LTE Co abundances from the current work and other 
source s. Closer to our NLTE results is the model of [Francois et al.l 
( I2OO4 ') with unadjusted metallicity-independent SN yields, but the 
offset between theoretical calculations and the spectroscopic abun- 
dances is still too large, ~ 0.25 dex. The chemical evolution 
models with me tallicity-depend ent SN yields jTinunes et al.lll993 : 
[Samland 1998; |g oswami & Pr antzos 200^ predict radically dif- 
ferent [Co/Fe] trends, which are ruled out both by NLTE and LTE 
sp ectroscopic ab undances of Co. The same is true for the results 
of lAmetll jlgTlh Co evolution calculations, whether with or with- 
out hydrostatic carbon bur ning. Th is discrepancy was also noted 
by [Xrnett et al. ( 197^ and lArnetj fl971), who proposed that the 
[Co/Fe] ratio should be strongly correlated with the SN la progen- 
itor metal content and the amount of mixing in the ISM. However, 
the effect that this would have on the abundances in metal-poor 
stars is unclear. 

The conclusion from the NLTE trend of Co in metal-poor stars 
is that massive stars overproduce Co relative to Fe. The observed 
decline of [Co/Fe] at [Fe/H] 1 can be explained by an underpro- 
duction of Co in SN la, but the possibility that SN II and SN la Co 
yields are metallicity-dependent is also not excluded. However, we 



6 CONCLUSIONS 

Deviations from LTE for Co i are controlled by overionization in 
the whole range of stellar parameters considered in this work. Thus, 
ionization equilibrium of Co is shifted towards lower number den- 
sities of Co I atoms compared to LTE. Excitation equilibrium of 
Co II is also disturbed due to optical line pumping. 

NLTE effects on the line formation occur for lines of both ion- 
ization stages. For Co i, reduced line opacity stipulated by the ove- 
rionization leads to general weakening of lines. Lines of Co 11 are 
affected by the deviation of line source functions from the Planck 
function. In the solar model atmosphere, these deviations are small, 
and it is safe to adopt LTE for Co 11. In the models with very 
low metallicities and supersolar temperatures, NLTE corrections 
for Co II lines are as large as -1-0.3 dex. NLTE abundance correc- 
tions for the lines of Co i vary from -l-O.l to -1-O.6 dex depending 
on the metallicity and the effective temperature. Our results for the 
giant models with low T^jf, low logg, and low [Fe/H] suggest that 
in atmospheres of giants deviations from LTE are larger than in 
dwarfs, although the main stellar parameter that controls the mag- 
nitude of NLTE effects in Co is metallicity. 

The solar abundance of Co derived from the NLTE analysis 
of Co I lines is 4.95 ± 0.04 dex, which is 0.13 dex higher than 
the LTE abundance. The abundance calculated from the single rel- 
atively unblended Co 11 line in the solar spectrum is 4.86 dex. 
The discrepancy of 0.09 dex between lines of two ionization stages 
points either at a failure of the NLTE approach to describe the ion- 
ization equilibrium of Co, or to an erroneous theoretical oscillator 
strength for the Co 11 transition. The NLTE abundances calculated 
with different scaling factors to cross-sections for H i collisions, 
0.05 < 5h < 5, are in agreement with the Co abundance measured 
in C I meteorites, 4.89 ± 0.01 dex. 

For the solar-metallicity stars, the incorrect treatment of hy- 
perfine structure leads to large errors in abundances. For the lines 
with saturated cores, the neglect of HFS leads to an abundance 
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overestimate by > 0.5 dex. Even the weak lines in the solar 
spectrum are not free from HPS effects, showing differences of 
0.05 . . .0.15 dex. The influence of HPS on the line profiles and 
abundances is negligible in the metal-poor stars; not exceeding 0.05 
dex. 

The difi^erential abundance analysis of 17 stars with -2.5 < 
[Pe/H] < demonstrated that the evolution of Co abundances in 
the Galaxy is radically different from all previously reported trends. 
The NLTE abundance ratios [Co/Pe] are supersolar for the thick 
disk stars, and at low metallicities [Co/Pe] are as large as 0.8 dex. 
This result suggests that Co is overproduced relative to Pe in short- 
lived massive stars. The GCE models, which utilize metallicity- 
dependent yields for SNe, predict a qualitatively different [Co/Pe] 
trend for the metal-poor stars, characterized by a plateau or [Co/Pe] 
depletion. We relate the discrepancy to the deficiency in SN yields. 
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Figure Al. Grotrian diagram of the Co i/Co ii model atom. Solid lines represent allowed and forbidden transitions included in the model atom. Ionization 
energies of Co i and Co n are 7.86 eV and 17. 1 eV, respectively. 
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Table Al. HPS constants A and B (in units of 10"^ cm"') representing magnetic dipole and electric quadrupole interactions for Co i and Co n levels. Level 
energies E are given in eV. The interaction constants for all Co i levels are taken from Pickering (1996). For the Co n levels, the data measured in this work 
are used. 
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